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We propose a simple encoding of charge-based quantum dot qubits which protects against fluctu- 
ating electric fields by charge symmetry of the logical states. We analyze the reduction of coupling 
to noise due to nearby charge traps and present single qubit gates. The relative advantage of the 
encoding increases with lower charge trap density. 
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I. INTRODUCTION 

Quantum computation faces considerable hurdles, one 
of the most serious being engineering physical systems 
performing coherent operations without the deleterious 
effects of decoherence^, particularly in the solid state. 
However by isolation and manipulation of states of quan- 
tum dot (QD) structures^, it may be possible to per- 
form many unitary operations within the dephasing time, 
a pre-requisite for quantum error correction (QEC) by 
means of Calderbank-Shor-Steane codesi 

Underlying logical QEC, a complementary strategy is 
to use Hilbert subspaces which couple least to noise pro- 
cesses, decoherence free subspaces (DFS^^ii. Practi- 
cal quantum computing will undoubtedly use elements of 
both. Charge-based QD quantum computing 8 &i2iii*i2ii^ 
is a prime candidate for DFS encoding as electric field 
coupling is a major source of decoherenc o 14 : 15 . Here, we 
present an architecture incorporating charge symmetry of 
the logical states to protect against electromagnetic fluc- 
tuations, analyze its resistance to charge trap noise and 
present single-qubit gates. Coupling to charge trap noise 
and decoherence is suppressed by several orders of mag- 
nitude compared to a conventional charge qubit, depend- 
ing on charge trap density. Alternatives to the passive 
control implied by DFS encoding include active control 
sequences, such as Bang-Bang controii&ii. 

In a typical charge-based QD qubit (Fig. [IJi) the posi- 
tion of an excess electron defines the logical states. Ide- 
ally, the logical states of the system should be eigenstates 
of the system Hamiltonian when the system is idle, i.e. 
tunnelling should be suppressed on practical time-scales 
by Vq. Furthermore, we assume that the system can be 
tuned, via Vo,i, such that the logical states are degen- 
erate, hence (known) relative dynamical phases can be 
neglected. 



II. CLASSICAL NOISE 

Fluctuations of the electromagnetic environment su- 
perimpose inhomogeneities on the potential seen by the 
charge states. An electric field component along the axis 




FIG. 1: a) 2QD Dipole Charge Qubit. Logical states, 
{|0}, |1}}, are defined by an excess electron on the left or right 
dot respectively. Symmetry electrodes (Vo,i) control a z ro- 
tations, and a barrier electrode (Vh) controls a x . An SET 
(omitted) measures in the logical basis, b) Temporally and 
spatially varying potential due to external field, c) 2-Electron 
4QD Quadrupole Qubit with |0) = a^jvac), |1) = fc T d T |vac). 

of the qubit will cause a sloping potential (Fig. ^p), in- 
ducing for each state a different dynamic phase, 

\m)=e-%ti^ dt '\j(0)), j =0,1 (1) 

where tj is the on-site energy fluctuation and q the elec- 
tron charge. Fluctuations drive superpositions \ip) = 
a|0) + f3\l) to mixed states, \ip){i>\ i-> |a| 2 |0)(0| + 
|/3| 2 |1)(1|. Furthermore, electrodes operating on nearby 
qubits will look like noise, i.e. it may only be practical 
to actively compensate for operations on nearest neigh- 
bors, but not those further away which may also cause 
unwanted perturbations. 

We generalize Eq. JIJ to multi-electron configurations, 
encoding logical states in many-particle states whose ge- 
ometry protects against decoherence (Fig. ^) 18 ' 19 . Two 
excess electrons in diagonally opposite dots define the 
logical states. Single square QDs in the limit of large dot 
size should display similar dynamics 20 21 . The 4QD ar- 
rangement has also been considered for Coherent Quan- 
tum Cellular Automata^ 2 * 2 ^, and for scalable qubitsiS. 
Measurement in the logical basis can be achieved by a 
single electron transistor (SET) adjacent to one of the 
dots in each qubiti^, or by using multiple SETs in a cor- 
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related mode 2 ^. Qutrits or higher dimensional systems 
may also be considered, e.g. a qutrit encoded as two elec- 
trons in a three dimensional 6QD octohedral structure. 

An external electric field induces phase shifts, as in 
Eq. QJ, where the energies to first order are 

eo(t) = e A + e c = 2y + ko-E+(-ko)-E = 2F 
d(t) = e B + e jD = 2V r + k 1 -E+(-k 1 )-E = 2t/, 

where V is the potential at the common centroid and 
£a,b,c,d are the on-site energy fluctuations of the re- 
spective QDs. The symmetrical distributions of charge 
ensure that each logical state acquires the same dynamic 
phase due to the external potential gradient. Thus, an 
initial superposition acquires an overall dynamic phase 
which is unobservable. 



III. CHARGE TRAP NOISE 

Though linear spatially varying potentials have no de- 
phasing effect on the 4QD qubit, charge trap fluctu- 
ators 2 ^ may pose a problemS5iSSiSLSSi2£. An occupied 
charge trap has a ~ 1/r potential, which perturbs the 
degeneracy of the DFS states. In principle the charge 
trap density can be made arbitrarily low but a few charge 
traps may be unavoidable in practice^, and charge trap 
noise may be a significant source of decoherence. 

To understand the effect of charge trap noise, consider 
a single charge trap coupled to the qubit via the Hamil- 
tonian H = k£(t)a z /2, where £ = ±1 is a Poisson process 
of rate A, and k is the coupling22i2i. Averaging over noise 
processes leads to a decay of the coherence of the qubit 
density operator, 



(poi(t))t = /MO) (e- ifc /o ^ 
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where lu = \/k 2 — A 2 . For many independent fluctuators 
with different rates Xj and couplings kj, the coherence 
decays in a non-Markovian manner (Fig. |2J, 



(poi(t)) e =poi(0)e-*£^ J] 
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A Taylor expansion of the solution J3J) about t — shows 
that the initial decay is parabolic 

<Poi)«/poi(0)«l-i 2 /2^fc, 2 + 0(t 3 ), t«l, (3) 



independent of Xj, and depends only on the effective 
coupling of the encoded qubits to the charge traps, 
k 2 cS = k 2 (Fig. Ok). Therefore, the short-term behav- 
ior will be dominated by the fluctuator that couples most 
strongly to the qubit, while the others mainly dampen 




FIG. 2: Qubit Coherence Decay. Initial state (|0) + |1))/V2 
coupled to 100 charge traps, = 12.46 x 10 9 h/s and 

feS* = 0.60 x 10 9 ?i/s. A total of 200 quantum trajectories 
were simulated and averaged. For all simulations, the 2QD 
qubit was 20nm long and the 4QD qubit was a 20nm square 
(e.g. P donors in Si, see Ref.H). Both were located 20nm be- 
low the layer in which the charge traps were located. Charge 
trap transition rate was 2 x 10 8 Hz. 



further oscillations of the coherence vector (Fig. |2J • Fur- 
thermore, the time it takes for the coherence to decay 
from 1 (maximal coherence) to p for p close to 1, which 
is of crucial importantance in quantum information pro- 
cessing, is inversely proportional to the effective coupling 
strength k eS , t p = ^2(1 - p)k~ 



c e g, and we have 



_(2) 



(2) 
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.(4)' 
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where k^ and k^ is the effective coupling strength for 
the two-dot and four-dot encoding respectively (Fig. |3| 
b). Thus, the ratio of the effective coupling strengths is 
a good measure for the superiority of the 4QD encod- 
ing versus the 2QD encoding — the former will be better 

provided that k^g /fe^i > 1, and the larger the ratio the 
greater the improvement. 

The 2QD and 4QD qubits couple differently to charge 



(2) traps. 



; (2) 



and 



kf OC Tj 



respectively, where 



is the distance between the qubit and each charge 
trap. The 4QD qubit has thus effectively a smaller "hori- 
zon" than the 2QD qubit. Hence, generally a charge trap 
would have to be situated closer to the 4QD qubit than 
a 2QD qubit to induce the same decoherence. Since the 
noise on the qubit is generally dominated by the clos- 
est fluctuator, whose typical distance is inversely propor- 
tional to defect density, the average relative effectiveness 
of the encoding is therefore expected to increase with de- 
creasing charge trap density, which is confirmed by com- 
puter simultations (Fig. HJ . 
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FIG. 3: Decoherence of 2QD and 4QD qubits. (a) The de- 
coherence decay times are inversely related to the effective 
coupling to fluctuators and show the same dependence for 
the 2QD and 4QD qubits. (b) The ratio of the short-term co- 
herence times for the 2QD and 4QD encodings are inversely 
proportional to the ratio of the effective coupling constants. 
Each point represents the average of 200 quantum trajectories 
of a qubit coupled to 100 randomly distributed fluctuators. 
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median of noise suppression 
— ■ linear fit p(x) = -9.412 x+ 100.646, 
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FIG. 4: Coupling and Noise vs Charge Trap Density. Noise 
suppression of the 4QD vs 2QD qubits was calculated for 
100 random charge trap distributions per density, which was 
then averaged. Coupling strength is in units of W 9 h/s, bars 
indicate the 10% — 90% ranges. At high densities, the mean 
charge trap spacing is comparable to the size of the qubit, 
leading to saturation effects. 



dipole is comparable to, and linear in the displacement, 
and given that fabricational precision should be at least 
a fraction of QD spacing for QIP purposes^, the extra 
dipole for the 4QD qubit should be much smaller than 
for a 2QD qubit, hence the encoding should still offer a 
noticeable advantage. 

To quantify the effect of asymmetry due to placement 
errors in the quantum dots, we performed extensive sim- 
ulations computing the effective couplings for various 
randomly perturbed 2QD and 4QD architectures for dif- 
ferent charge trap densities and a wide range of charge 
distributions. The simulations show that for reasonable 
errors (~ 10% placement error), the efficiency of the 
scheme is only modestly affected over a wide range of 
fluctuator densities (Fig. J5J . 



IV. ROBUSTNESS OF ENCODING 



The results in the previous section show that the 4QD 
encoding can substantially increase short-term coherence 
times for an ideal geometry. However, any physical im- 
plementation is likely to deviate from the perfect sym- 
metry of the ideal quantum dot structure. The scheme's 
sensitivity to such deviations is thus an important prac- 
tical consideration. 

Non-ideal geometry, e.g. due to imperfect QD place- 
ment, will introduce a dipole moment, spoiling decou- 
pling from external fields and reducing robustness to 
charge trap noise. However, as the magnitude of this 



V. QUANTUM GATES 

We now consider implementing a universal set of quan- 
tum gates, {of,(Tx ,c — erj}. Ideally, we would like 
all states involved during gate operations to belong to 
the DFS. This suggests an adiabatic holonomic control 
scheme^*^^. However, the requirement of additional 
quantum dots for generating holonomies, charge symme- 
try constraints on auxiliary dot positioning, and the com- 
plexity of pulse sequences all offset possible advantages 
of holonomic control. Alternatively, rapidly modulating 
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FIG. 5: Decoupling vs Placement Error. We simultaneously 
and independently purturbed the positions of all quantum 
dots by a Gaussian displacement with standard deviation a 
(expressed as a fraction of the array side-length). This was 
repeated 1000 times for each of 500 fluctuator distributions. 
We plot of median w.r.t fluctuator distributions of the mean 
over the perturbations of the decoupling ratios as a function 
of error magnitude for different effective fluctuator densities. 
Even with a 10% displacement error (i.e. a — 0.1), the 4QD 
qubit is still effective. 



tunneling between dots can implement the required log- 
ical gates quickly. If the gate time is short enough, tran- 
sient population in non-DFS states, should have minimal 
coupling to electric field fluctuations. Intra-dot charging 
should effectively suppress double occupation, which can 
be further enhanced by ensuring all spins are parallel so 
that each orbital cannot have more than one electron. 

We describe the four dot system as a two-electron, 
four-site Hubbard model. The electronic creation oper- 
ators are defined a\ tf, c<, eft for dots A, B, C, D respec- 
tively, labelled in clockwise fashion (Fig-QJ;). Firstly, the 
phase gate (a rotation of the Bloch sphere by angle 4> 
around the z-axis) of = diag(l, e 1 ^) is achieved by bias- 
ing one pair of diagonally opposite quantum dots with 
respect to the other, 



m* [Vow-v^df, 



where Vq , V\ are the on-site potentials of the quantum 
dots defining the |0) = a^ct|vac) and |1) = weft\va,c) 
states respectively. 

Next, the a x gate requires inter-dot tunnelling. We 
allow tunnelling between dots A <-> D and B <-> C. Al- 
lowing A <-> B and C <-> D tunneling as well leads to 
similar dynamics but at the expense of extra control elec- 
trodes and more non-DFS states involved. With vertical 
tunneling only, the available state space is spanned by 
four states, {|0), |eo) = a^tf |vac), |ei) = c^d^vac)}. 
Hence the Hamiltonian with no tunneling is, Hq = 
diag(0, 0, S, 5) in the above basis, and where we have 
taken the (degenerate) ground state energy to be and S 
is the energy of the non-diagonal states {|eo), |ei)} due to 



Coulomb repulsion of the two electrons. We now switch 
on equal tunneling with rate O in the vertical direction, 



H, 



tunnel 
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For convenience, we normalize S 
to this. The eigenstates of H to t 



■ 1, and scale f2 relative 
#o + -^tunnel are 
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where the eigen-energies are E\ 

Vi + im 2 )/2,E A ~ 



Q,E 2 



(1 + VI + 16ft 2 )/2, E 4 = (l- Vl + 16f2 2 )/2. Tunnelling 
between dots mixes the states so that {|0), |1)} are no 
longer eigenstates of H to t- Transitions between |0) and 
|1) cannot occur directly but only via transient occupa- 
tion of the non-DFS states. 

In order to achieve a n/2 rotation around an axis 
(cos 7, sin 7, 0) lying on the equator of the Bloch sphere 

(which is equiavalent to a a^ 2 gate up to rotations) 
|0> 1 ^ |0) + e*T|l), |1) 1 — ^ |0) - e^|l), we require E 3 
and E4 to be ratio nal. This leads to the conditions, 
4f2 = y/ (n/m) 2 — 1 where {n/2,m} C Z + , n > m and 
gcd(n, m) — 1. These requirements derive from the fact 
that the amplitudes of the |0) and |1) should be equal in 
magnitude when the amplitudes of the non-DFS states 
are zero, leading to jm/n = 1/2 + k, {j, k} C Z + , and 
the gate time tt = 2jrmt/n. When tj = 7rm, we achieve 
the operation with 7 = ir(n — m)/2. If tf — 2irm, we 
perform a logical NOT (|0) <-> |1». 

The minimum gate time for a 7r/2-gate is tf — tt when 
2 < n even and m = 1 (Fig. [fjjl. If m > 1, the time re- 
quired to implement the gate increases. Coulomb repul- 
sion favors the diagonal charge configurations but tran- 
sient population in the other states will still occur. For 
n/m — > 1, the gate time is on the order of 7rm but the 
maximum transient population scales as (n 2 — m 2 )/n 2 . 
The integrated population in the non-DFS states during 
the total gate time is proportional to m(n 2 — m 2 )/n 2 and 
thus using smaller tunneling rates does not improve the 
overall transient occupation of non-DFS states. 

The average gate error, £ = 1 — T where T is the 
average fidelitySi*2I, for different charge trap couplings 
(densities) was simulated for 4QD and 2QD qubits. The 
ratio of the errors, presented in Fig. [71 show that despite 
transient population in non-DFS states during the oper- 
ation of the a x -gate, the 4QD configuration still shows 
a significant advantage over the 2QD qubit. 

A universal two-qubit controlled-phase (c— <p) gate may 
be implemented as suggested in earlier work2&. A tran- 
sient deformation of the charge distribution of adjacent 
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qubits by the use of auxiliary quantum dots would allow 
modulation of an effective a z ® a z interaction. Charge 
symmetry could still be maintained during the gate by 
use of auxiliary dots. 



VI. CONCLUSION 
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FIG. 6: Populations During cr^^-Gate vertical tunneling. 
An initial state |0) is transformed into (|0) + i|l})/\/2 with 
transient population in non-DFS subspace |e). a) fi large: 
n = 2, to = 1, b) f2 small: n = 20, m = 19. 
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2 — 1. Fidelity was calculated from 50 trajectories 
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per initial state { | ± a;} , j ± y) , j ± z) } 



We have analyzed the noise suppression of a 2-electron 
4QD qubit encoding, which decouples from linearly vary- 
ing fields. For nearby charge fiuctuators, the decoupling 
depends upon the exact distribution, but analytic and 
numerical results show considerable enhancement of noise 
resistance and coherence times, increasing at low charge 
trap densities. This advantage persists in the presence 
dot placement errors. We show how to construct sin- 
gle qubit rotations, and two-qubit gates may be pos- 
sible via previously proposed schemes for conventional 
charge qubits. Characterization and tuning of the 4QD 
qubit should be be similar to that for a conventional 2QD 
qubit. The architecture requires only a modest increase 
in complexity and may also be applied to systems such 
as super-conducting charge qubits. 
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